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Abstract 

Measurements of dijet azimuthal decorrelations in pp collisions at v^s = 7 TeV us- 
ing the CMS detector at the CERN LHC are presented. The analysis is based on an 
inclusive dijet event sample corresponding to an integrated luminosity of 2.9 pb^^. 
The results are compared to predictions from perturbative QCD calculations and var- 
ious Monte Carlo event generators. The dijet azimuthal distributions are found to be 
sensitive to initial-state gluon radiation. 
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High-energy proton-proton collisions with high momentum transfer are described within the 
framework of Quantum Chromodynamics (QCD) as point-like scatterings between the proton 
constituents, collectively referred to as partons. The outgoing partons manifest themselves, 
through quark and gluon soft radiation and hadronization processes, as localized streams 
of particles, identified as jets. At Born level, dijets are produced with equal transverse mo- 
menta Pt with respect to the beam axis and back-to-back in the azimuthal angle (Acp^ijet = 
l^'jeti — 9]et2\ = tt). Soft-gluon emission will decorrelate the two highest pj (leading) jets and 
cause small deviations from n. Larger decorrelations from n occur in the case of hard multijet 
production. Three-jet topologies dominate the region of 27r/3 < A^dijet < tz, whereas angles 
smaller than In/ 3 are populated by four-jet events. 

Dijet azimuthal decorrelations, i.e., the deviation of Aipdijet from tt for the two leading jets 
in hard-scattering events can be used to study QCD radiation effects over a wide range of 
jet multiplicities without the need to measure all the additional jets. Such studies are impor- 
tant because an accurate description of multiple-parton radiation is still lacking in perturbative 
QCD (pQCD). Experiments therefore rely on Monte Carlo (MC) event generators to take these 
higher-order processes into account in searches for new physics and for a wide variety of preci- 
sion measurements. The observable chosen to study the radiation effects is the differential dijet 
cross section in A^dijet/ normalized by the dijet cross section integrated over the entire A^dijet 
phase space, (l/c^dijet) (^c^dijet/ ^^A^dijet)- By normalizing the A^dijet distributions in this manner, 
many experimental and theoretical uncertainties are significantly reduced. Measurements of 
dijet azimuthal decorrelations at the Tevatron have previously been reported by the DO collab- 
oration Ijll. In this Letter, we present the first measurements of dijet azimuthal decorrelations 
in pp collisions at a/s = 7 TeV at the CERN Large Hadron Collider (LHC). 

The central feature of the Compact Muon Solenoid (CMS) apparatus is a superconducting 
solenoid, of 6 m internal diameter, providing an axial field of 3.8 T. Charged particle trajec- 
tories are measured by the silicon pixel and strip tracker, covering < < 2 tt in azimuth 
and \r]\ < 2.5, where pseudorapidity rj = — In [tan(0/2)] and 9 is the polar angle relative 
to the counterclockwise proton beam direction with respect to the center of the detector. A 
lead-tungstate crystal electromagnetic calorimeter and a brass/ scintillator hadronic calorime- 
ter surround the tracking volume. The calorimeter cells are grouped in projective towers of 
granularity At] x Aq) = 0.087 x 0.087 at central pseudorapidities. The granularity becomes 
coarser at forward pseudorapidities. A preshower detector made of silicon sensor planes and 
lead absorbers is installed in front of the electromagnetic calorimeter at 1.653 < < 2.6. 
Muons are measured in gas-ionization detectors embedded in the steel magnetic field return 
yoke. A detailed description of the CMS detector can be found elsewhere [2J. 

CMS uses a two-tiered trigger system to select events online; Level-1 (LI) and the High Level 
Trigger (HLT). In this analysis, events were selected using two inclusive single-jet triggers that 
required an LI jet with > 20 GeV (30 GeV) and an HLT jet with pr > 30 GeV (50 GeV). 
The jets at LI and HLT are reconstructed using energies measured by the electromagnetic and 
hadronic calorimeters and are not corrected for the jet energy response of the calorimeters. The 
trigger efficiency for a given corrected pr threshold of the leading jet (pr'"^") was measured 
using events selected by a lower-threshold trigger. For the event selection, pj"^^^ thresholds 
were chosen so that this efficiency exceeded 99%. The corresponding offline corrected pf""^" 
values are 80 GeV (110 GeV) for the low (high) threshold jet trigger. 

Jets were reconstructed offline using the anti-Zcj clustering algorithm with a distance parame- 
ter R = 0.5 The four-vectors of particles reconstructed by the CMS particle-flow algorithm 
were used as input to the jet-clustering algorithm. The particle-flow algorithm combines in- 
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formation from all CMS sub-detectors to provide a complete list of long-lived particles in the 
event. Muons, electrons, photons, and charged and neutral hadrons are reconstructed indi- 
vidually. As a result, the residual corrections to the jet four-vectors, arising from the detector 
response, are relatively small (at the level of 5-10% in the central region) |^. A detailed de- 
scription of the particle-flow algorithm can be found elsewhere [51161. 

Spurious jets from noise and non-collision backgrounds were eliminated by applying loose 
quality cuts on the jet properties |7J. Events were required to have a primary vertex recon- 
structed along the beam axis and within 24 cm of the detector center |8|. Further cuts were 
applied to reject interactions from the beam halo. Events were selected having two leading jets 
each with pj > 30 GeV and rapidity \y\ < 1.1, where y = jln [(E + p^) / (E — p^)], with E 
being the total jet energy and pz the projection of the jet momentum along the beam axis. Each 
event is put into one of five mutually exclusive regions, which are based on the pj"^^^ in the 
event. The five regions are: 80 < pr"''"' < 110 GeV, 110 < pr'"'"' < 140 GeV, 140 < pr'"'"' < 200 
GeV, 200 < pr""'"' < 300 GeV, and 300 GeV < pj"'^''. The data correspond to an integrated lu- 
minosity of 0.3 pb^^ for the lowest pj^^^ region and 2.9 pb^^ for the other p^"^^" regions. The 
uncertainty on the integrated luminosity is estimated to be 11% |9|. After the application of all 
selection criteria, the numbers of events remaining in each of the five p^"^^" regions, starting 
from the lowest, are: 60837, 160388, 69009, 14383, and 2284. 

The A^dijet distributions are corrected for event migration effects due to the finite jet pj and 
position resolutions of the detector. The distributions are sensitive to the jet pr resolution be- 
cause fluctuations in the jet response can cause low-energy jets to be misidentified as leading 
jets, and events can migrate between different pj"^^^ regions. The finite resolution in azimuthal 
angle causes event migration between Acp^ijet bins, while the resolution in rapidity can move 
jets in and out of the central rapidity region {\y\ < 1.1). The correction factors were determined 
using two independent MC samples: pythia 6.422 (PYTHIA6) [lOJ tune D6T illj|, and her- 
WIG++ 2.4.2 [12]. The pj, rapidity, and azimuthal angle of each generated jet were smeared 
according to the measured resolutions [131 . The ratio of the two dijet azimuthal distributions 
(the generated distribution and the smeared one) determined the unfolding correction factors 
for each pj"^^^ region, for a given MC sample. The average of the correction factors for each 
p^max region from the two MC samples were used as the final unfolding correction applied to 
data. The unfolding correction factors modify the measured A^dijet distributions by less than 
2% for 57r/6 < A^^ijet < ^- For Aipdijet ~ 7r/2, the changes range from —11%, for the highest 
p^max region, to —19%, for the lowest. 

The main sources of systematic uncertainty arise from uncertainties in the jet energy calibra- 
tion, the jet Pt resolution, and the unfolding correction. The jet energy calibration uncertainties 
have been tabulated for the considered phase space in the variables of jet pr and rj [4J. Typ- 
ical values are between 2.5% and 3.5%. The resulting uncertainties on the normalized A^dijet 
distributions range from 5% at A^dijet ~ 7r/2 to 1% at A(pdijet ~ ^- The effect of jet pj reso- 
lution uncertainty on the A^dijet distributions was estimated by varying the jet pj resolutions 
by ±10% Ills ! and comparing the A^dijet unfolding correction before and after the change. This 
yields a variation on the normalized A^dijet distributions ranging from 5% at A^dijet ~ 7r/2 
to 1% at A^dijet ~ ^- The uncertainties on the unfolding correction factors were estimated by 
comparing the corrections from different event generators and PYTHIA6 tunes that vary sig- 
nificantly in their modelling of the jet kinematic distributions and A^dijet distributions. The 
resulting uncertainty varies from 8% at A^dijet ~ 7r/2 to 1.5% at A(pdijet ~ tt. The systematic 
uncertainty from using a parametrized model to simulate the finite jet pj and position resolu- 
tions of the detector to determine the unfolding correction factors was estimated to be about 
2.5% in all p;^™^" regions. The combined systematic uncertainty, calculated as the quadratic 
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sum of all systematic uncertainties, varies from 11% at A^^ijet ~ 7r/2 to 3% at A^dijet ~ tt. 

The corrected differential A^dijet distributions, normalized to the integrated dijet cross section, 
are shown in Fig. [ijfor the five p^"^^" regions. The distributions are scaled by multiplicative 
factors for presentation purposes. Each data point is plotted at the abscissa value for which the 
predicted differential A^dijet distribution has the same value as the bin average obtained using 
PYTHIA6 tune D6T, which provides a good description of the data ||T4| . 

The A^dijet distributions are strongly peaked at n and become steeper with increasing pj'"^". 
The simulated A^dijet distributions from the PYTHIA6 (D6T and Z2 llSl tunes), pythia 8.135 
(PYTHIA8) llll, HERWIG++, and MADGRAPH 4.4.32 [17] event generators are presented for 
comparison. The MADGRAPH generator is based on leading-order matrix element multiparton 
final-state predictions, using PYTHIA6 for parton showering and hadronization, and the MLM 
method ||T8| to map the parton-level event into a parton shower history. The MADGRAPH pre- 
dictions included tree-level processes of up to four partons. For PYTHIA6, pythiaS, and MAD- 
GRAPH event generators the CTEQ6L [|19J parton distribution functions (PDFs) were used; for 
HERWIG++, the MRST2001 PDFs Il20l . 

Figure |2] shows the ratios of the measured A^dijet distributions to the predictions of PYTHIA6, 
PYTHIA8, HERWIG++, and MADGRAPH in the five pj^^^ regions. The combined systematic un- 
certainty on the experimental measurements is shown by the shaded band. The predictions 
from PYTHIA6 and HERWIG++ describe the shape of the data distributions well, while MAD- 
GRAPH (pythiaS) predicts less (more) azimuthal decorrelation than is observed in the data. 

Figure |3| displays a comparison between the measured A^dijet distributions and the predic- 
tions ofpQCD calculations from the parton-level generator NLOJET++ [21J within the FASTNLO 
framework ||22| . The predictions near Acpdijet = have been excluded because of their sensi- 
tivity to higher-order corrections not included in the present calculations. The leading-order 
(LO) curves represent processes with three partons in the final state, normalized to the LO cTdijet 
(2 2) cross section. The next-to-leading-order (NLO) predictions include 2 — > 3 processes at 
NLO, normalized to (T^i]ei at NLO: 



C^'dijet 



dcTdijet 

X 

(N)LO rfAcpdijet 



(N)LO 



Uncertainties due to the renormalization {}ir) and factorization {ji^) scales are evaluated by 
varying the default choice of = }if = pr'"^" between pt"^^''/2 and 2pj^^^ in the following 
six combinations: {}lrr^f) = (Pt"'^V2, Pt°'^V2), {Ipj"'^'', Ipr"^^^), (pT"'^^ Pt"'^V2), (pT'"''^ 
2p^max^^ |-p^maxy2^ p^max-j^ ^2pj,max^ p^max-j 'pj-jggg gcale Variations modify the predictions 

of the normalized Acpdijet distributions by less than 50%. The PDFs and the associated uncer- 
tainties were obtained from CTEQ6.6 [19J. The PDF uncertainties were derived using the 22 
CTEQ6.6 uncertainty eigenvectors and found to be 9% at A^dijet ~ ^/2 and 2% at A^dijet < ^• 
Following the proposal of the PDF4LHC working group [23 1, the impact of other global PDF 
fits [|24ti26l was investigated and found to be negligible in the context of this analysis. 

Non-perturbative corrections due to hadronization and multiple-parton interactions were ap- 
plied to the pQCD predictions. The correction factors were determined from the PYTHIA6 
and HERWIG++ simulations and modify the predictions from +4% (Acpdijet ~ ^r) to —13% 
(A^dijet ~ 7t/2). The uncertainty due to the non-perturbative corrections is estimated to be 
6% at A^dijet ~ 7^/2 and 2% at Ag?dijet ~ /r.The ratios of the measured Ag?dijet distributions to 
the NLO pQCD predictions are shown in Fig. |4] The effect due to the scale variations, as well 
as the uncertainties due to PDFs and non-perturbative corrections, are also shown. The NLO 
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predictions provide a good description of the shape of the data distributions over much of the 
A^'dijet range. Compared to the data, the reduced decorrelation in the theoretical prediction and 
the increased sensitivity to the fir and }if scale variations for A^^ijet < 2n/3 shown in Fig. |4] 
are attributed to the fact that the pQCD prediction in this region is effectively available only at 
leading order, since the contribution from tree-level four-parton final states dominates. 

The sensitivity of the A^^ijet distributions to initial-state parton shower radiation (ISR) is in- 
vestigated by varying the input parameter /cjgR (PARP(67)) in PYTHIA6 tune D6T. The product 
of kisR and the square of the hard-scattering scale gives the maximum allowed parton virtu- 
ality (i.e., the maximum allowed pj) in the initial-state shower. Previous studies have shown 
that fcisR is the only parameter in PYTHIA6 that has significant impact on the A^^ijet distribu- 
tions llZl. The default value of kisR in PYTHIA6 tune D6T is 2.5, determined from the DO dijet 
azimuthal decorrelation results [IJ. Figure |5] shows comparisons of the measured Acpjijet dis- 
tributions to PYTHIA6 distributions with various kisR values. The effects are more pronounced 
for smaller Acp^^et angles, where multi-gluon radiation dominates. Varying kjsR by ±0.5 about 
its default value yields a change of about 30% on the PYTHIA6 prediction for A^^ijet ~ 7r/2, 
suggesting that our results could be used to tune parameters in the MC event generators that 
control radiative effects in the initial state. In PYTHIA6 tune D6T, the maximum pj allowed in 
final-state radiation parton shower is controlled through the parameter parp(71). We varied 
the value of this parameter from 2.5 to 8 (the default value is 4.0) and observed less than ~ 10% 
changes in the A^dijet distributions in all pj regions. 

In summary, we have measured dijet azimuthal decorrelations in different leading-jet pj re- 
gions from pp collisions at v^s = 7 TeV. The PYTHIA6 and HERWIG++ event generators are 
found to best describe the shape of the measured distributions over the entire range of A^dijet- 
The predictions from NLO pQCD are in reasonable agreement with the measured distributions, 
except at small Acpdijet where multiparton radiation effects dominate. The A^dijet distributions 
are found to be sensitive to initial-state gluon radiation. 
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Figure 1: Normalized A^^ijet distributions in several p^"^^" regions, scaled by the multiplica- 
tive factors given in the figure for easier presentation. The curves represent predictions from 
PYTHIA6, PYTHIAS, HERWIG++, and MADGRAPH. The error bars on the data points include 
statistical and systematic uncertainties. 
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Figure 2: Ratios of measured normalized A^dijet distributions to PYTHIA6, PYTHIA8, HER- 
WIG++, and MADGRAPH predictions in several pj"^^^ regions. The shaded bands indicate the 
total systematic uncertainty. 
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Figure 3: Normalized A^^ijet distributions in several p^"^^" regions, scaled by the multiplicative 
factors given in the figure for easier presentation. The curves represent predictions from LO 
(dotted line) and NLO pQCD (solid line). Non-perturbative corrections have been applied to 
the predictions. The error bars on the data points include statistical and systematic uncertain- 
ties. 
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Figure 4: Ratios of measured normalized A^^ijet distributions to NLO pQCD predictions with 
non-perturbative corrections in several pj"^^^ regions. The error bars on the data points in- 
clude statistical and systematic uncertainties. The effect on the NLO pQCD predictions due 
to ]!,■ and ji f scale variations and PDF uncertainties, as well as the uncertainties from the non- 
perturbative corrections are shown. 
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Figure 5: Ratios of measured normalized A^^ijet distributions to PYTHIA6 tune D6T with var- 
ious values of /cjsr in several pj^^^ regions. The shaded bands indicate the total systematic 
uncertainty. 
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